Free-volume kinetic models of granular matter 
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We show that the main dynamical features of granular media can be understood by means of 
simple models of fragile-glass forming liquid (Kob and Andersen, Phys. Rev. E 48 4364, 1993) 
provided that gravity alone is taken into account. In such lattice-gas models of cohesionless and 
frictionless particles, the compaction and segregation phenomena appear as purely non-equilibrium 
effects unrelated to the Boltzmann-Gibbs measure which in this case is trivial. They provide a 
natural framework in which slow relaxation phenomena in granular and glassy systems can be 
explained in terms of a common microscopic mechanism given by a free- volume kinetic constraint. 



Vibrated powders are one of the most interesting exam- 
ples of non-equilibrium systems where excluded-volumc 
interactions are known to play a crucial role. Based on 
the free-volume notion the analogy with glassy dynam- 
ics (and its limit) was early recognized JD. Indeed, the 
static and dynamic properties of powders seem to be even 
more intriguing than those appearing in amorphous sys- 
tems |^-^]. A natural question that arises is to what 
extent a powder can be considered as a glass . 

In this Letter we show that compaction and segrega- 
tion phenomena can be reproduced within a simple mi- 
croscopic model of fragile-glass forming liquid ||, pro- 
vided that gravity alone is taken into account. The 
key ingredient of the model is a free- volume kinetic con- 
straint involving only a selection of the possible configu- 
ration changes compatible with detailed balance and the 
Boltzmann-Gibbs distribution. Such kinetic constraints 
are so effective in describing slow relaxation phenomena 
that there is no need to invoke either a specific ener- 
getic interaction or a particular geometry of particles. 
Consequently, the ground state structure and the ther- 
modynamics of the model turn out to be trivial. This 
is the most peculiar feature of this lattice-gas model as 
compared to the previously proposed ones (see, for ex- 
ample JT^hTJ]). It follows that compaction and segre- 
gation phenomena are, at least in this case, purely dy- 
namical effects unrelated to the Boltzmann-Gibbs mea- 
sure. In particular, we find that the subtle interplay 
of kinetic constraints and the gravity driving force can 
account for several experimental findings such as loga- 
rithmic compaction, vibration dependence of the asymp- 
totic packing density and segregation phenomena Jl5| , ^6| . 
Such effect appear below a well defined vibration ampli- 
tude corresponding to the dynamical jamming transition 
of the model. Of course, the dynamics we consider is 
not intended to represent in a realistic way the micro- 



scopic granular dynamics but rather to provide an effec- 
tive stochastic description, typically at mesoscopic level. 

The model we consider was introduced by Kob and 
Andersen as a lattice-gas model of fragile-glass and used 
to test the predictions of the mode-coupling theory ||. 
In its original version the system consists of N particles 
on a simple cubic lattice where there can be at most one 
particle per site. There is no cohesion energy among par- 
ticles: the Hamiltonian is zero, 7i = 0. At each time step 
a particle and one of its neighbouring sites are chosen 
at random; the particle can move to the new site if it is 
empty and if the particle has less than v nearest neigh- 
bours before and after it has moved. This kinetic rule is 
time-reversible and therefore the detailed balance is sat- 
isfied. At high density the dynamics slows down because 
the reduced free-volume makes harder for a particle to 
meet the constraints. Above a certain threshold density 
the particles are so interlocked that macroscopic struc- 
tural rearrangements are no longer possible and the mo- 
bility steeply falls to zero. Such a mechanism is enough 
to provide the basic glassy phenomenology J|,[l7],[l8| • 

In the following we consider the Kob- Andersen model 
in presence of gravity. The system is described now by 
the Hamiltonian 



H = mg hi 



(1) 



where g is the gravity constant, rij = 0, 1 are the occupa- 
tion variables of the particles and hi is their height. We 
set throughout the mass of particles to one, m = 1. Parti- 
cles are confined in a box closed at both ends and with pe- 
riodic boundary condition in the horizontal direction. We 
then assume that the random diffusive motion of "grains" 
produced by the mechanical vibrations on the box can be 
modeled by a thermal bath at temperature T: particles 
satisfying the kinetic constraints can move according to 
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the Metropolis rule with probability min[l, x~ ], where 
Ah = ±1 is the vertical displacement in the attempted 
elementary move and x = exp(— g/T) represents the "vi- 
bration". We set the constraint threshold at v = 5. The 
system underlying geometry is that of a body centered 
cubic lattice. It gives the advantage of using a single pa- 
rameter to describe the combined effect of vibration and 
gravity. In this case the Markov process generated by 
the kinetic rules is irreducible on the full configurations 
space provided that the box height is large enough: in- 
deed it is always possible to find a path connecting any 
two configurations by letting, e.g., particles expand in the 
whole box. Therefore the static properties of the model 
are those of a non-interacting lattice-gas in a gravity field 
and they can be easily computed. We have investigated 
by extensive Montecarlo (MC) simulation the dynami- 
cal behavior of this model for a system of height 16L and 
transverse surface L 2 /2 with L = 20, for a number of par- 
ticles N — 16000. The runs are carried out over 10 6 MC 
time sweeps (MCs) and the observables averaged typi- 
cally over 20 realizations of dynamical evolution for each 
vibration amplitude x. The packing density is computed 
in the lower 25% of the system. 

When particles are subject to gravity and no vibra- 
tions, x — 0, the model has a single ground state with 
equilibrium packing density p = 1. It could be attained, 
e.g., by letting the particles fall from the top of the box 
Therefore we prepare the system in a 
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one at a time 

fluffy initial state by placing the particles at random in 
the upper half of the box and letting them fall randomly, 
x = 0, until they are not able to move anymore. We can 
assume that the mechanically stable state produced by 
this extensive statistical handling of particles represents 
a random loose packed state. In this state the average 
packing density turns out to be p r i p ~ 0.707. Once the 
system is prepared in a fluffy configuration the vibration 
with amplitude x is turned on. We used a continuous 
shaking procedure which is computationally more conve- 
nient and gives results qualitatively similar to the discon- 
tinuous tapping. In fig. [^a we plot the packing density, 
p{t) , as a function of time for several values of x in the 
full range of weak and strong tapping. In spite of the 
relatively high value of p r i p and the fact that kinetic con- 
straints can be thought to mimic an effective short-range 
repulsion, the packing density slowly increases as time 
goes on and its asymptotic value is an increasing func- 
tion of the vibration amplitude, even for quite high value 
of x, as in typical compaction experiments Jig] . The slow 
increase of p(t) is consistent, in the full range of weak and 
strong tapping, with a function of the form |fl5j: 



p(t) = p 



1 + B In (1 + t/r) ' 



(2) 



(see fig. [j]b). This inverse logarithmic law was first discov- 
ered in the experiments of Chicago group jlBI] , and is well 



supported by numerical simulations |12}-|l4|] , and analyt- 
ical approaches ]20|-p2]] essentially based on free-volume 
models. However, we find that in an intermediate range 
of vibration a power-law fit also gives good results, while 
for stronger vibration, x > 0.3, a modified logarithmic 
law works much better; a detailed discussion of these re- 
sults will be presented elsewhere. We find that the fit 
parameters, poo, B and r all depend on x. In particular, 
the asymptotic packing density, poo (x) , is always smaller 
than one and it depends smoothly on x in agreement 
with the experimental results p5[ |. There is an optimal 
vibration amplitude, x max ~ 0.9, at which p^ reaches 
a maximum; above cc max the asymptotic density become 
closer and closer to the equilibrium value (see fig. |l|c). 
The knowledge of thermodynamic properties allow us to 
find the point x g above which the dynamical results co- 
incide with the equilibrium ones. We find x s ~ 0.96 
(see fig. [|c). This critical vibration amplitude defines 
the jamming transition of the model which is the direct 
analogue of the glass transition. Below x g the system is 
no longer able to equilibrate on the experimentally ac- 
cessible time-scales. This is essentially due to a kinetic 
bottleneck at high packing density the number of paths 
leading to the equilibrium configurations is much smaller 
than the one leading elsewhere. Therefore, even if the 
dynamics is in principle always able to reach equilibrium 
configurations, actually for x < x g such configurations 
are not effectively accessible. We have checked that the 
usual MC procedure, obtained by removing the kinetic 
constraints, leads to the correct equilibrium configura- 
tion for any x. In this case the equilibrium state is easily 
reached and the corresponding packing density p e q{x) is 
a decreasing functions of x. Paradoxically, Poo(x) and 
p C q(x) exhibit opposite monotony properties as long as 

X <C ^max- 

We now turn to the behavior of the global density pro- 
file. In fig. H we can observe that near the bottom wall 
the strong packing effects give rise to highly structured 
density profiles, with exponentially damped oscillations 
over a length scale which decreases with the vibration 
strength. This effect is interestingly similar to the lay- 
ering phenomena appearing in simple fluids in confined 
geometry which are usually interpreted in terms of the 
"missing neighbours" interaction of the particles near 
the wall p3|j . In our case the particle- wall interaction is 
purely kinetic: the "particles" building up the wall can 
affect or not the number of nearest neighbours particles 
in the kinetic rule. We checked that in both cases the 
effect persists and results in a more or less dense bottom 
layer according to the "attractive" or "repulsive" nature 
of the kinetic particle-wall interaction. An experimental 
investigation of this point would be welcome. The expo- 
nential damping of these oscillations could explain why 
this fluid-like behavior is hardly observed in the experi- 
ments. A similar effect takes place near the top of the 
granular packing. However, due to the roughness of the 
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free surface the effect is less pronounced. 

We finally consider the behavior of a mixture with dif- 
ferent kinetic constraints. Indeed, size and mass segrega- 
tion phenomena are another puzzling behavior exhibited 
by powders and several distinct mechanisms have been 
proposed - 26 1 . Our main concern here is to show 

that, in their simplest form, segregation phenomena also 
appear in the present model. To keep things as simple 
as possible we consider only the case of a binary mix- 
ture with particles of identical mass. The only difference 
between the two components is purely kinetic: one type 
of particle is constrained while the other one is not; the 
latter particles can however block the former ones. The 
Hamiltonian of the system is given by (0), and the mi- 
croscopic reversibility is still satisfied. In hg. || we plot 
the time evolution of the density profile of a binary mix- 
ture starting from a random homogeneous distribution of 
particles. Although the packing fraction of both compo- 
nents increase, the system undergoes, at the same time, a 
spontaneous unmixing with the constrained particles seg- 
regating at the top of the mixture (see fig. ||). This turns 
out to be the case even when the mass of constrained 
particles is bigger than the unconstrained ones, at least 
in a certain range of the mass ratio of the two compo- 
nents. Such effect can be simply explained in terms of the 
different mobilities: during the shaking process the un- 
constrained "small" particles easily fill the gap beneath 
the constrained "large" particles whenever they open up; 
this decreases the probability of the slower constrained 
particles to fall down. A similar percolation or sifting 
mechanism was found in ref . p7[ . In our case it is clearly 
not related to an equilibrium phase transition, since the 
Hamiltonian is "blind" to the constraints. 

To summarize, we have investigated some aspects of 
granular dynamics in a three dimensional gravity-driven 
lattice-gas model. The key ingredient of the model is a 
free-volume constraint implemented in a purely kinetic 
way |§. Below a critical value of vibration the sys- 
tem undergoes a purely dynamical jamming transition. 
For weak vibrations the packing density logarithmically 
increases towards an asymptotic value which depends 
on the vibration amplitude, in agreement with experi- 
ments |jj"5|| . Near the bottom wall and the free surface 
of the granular cluster the strong packing effects induce 
highly structured density profiles; this phenomenon could 
be closely related to the one appearing in confined liquids. 
When a mixture of particles with different constraints is 
considered a vibration-induced phase-separation appear 
in the system. We emphasize that the results mentioned 
above are obtained in a model with trivial thermody- 
namics without introducing any form of quenched ran- 
domness neither in the energetic interaction nor in the 
shape of particles. Therefore all interesting effects have a 
purely dynamical nature and manifest themselves as a de- 
parture from the Boltzmann-Gibbs equilibrium measure. 
It would be interesting to derive these non-equilibrium 



features in terms of the Edwards' measure [E8[. 
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FIG. 1. Compaction experiment data: a) packing density 
p(t) computed in the lower 25% of the system vs time t for sev- 
eral values of the vibration amplitude x. The initial configura- 
tion is prepared in a random loose packed state (p r i p ~ 0.707); 
b) Logarithmic fit of compaction data in the weak (x < 0.3) 
and strong (a; > 0.3) tapping regime. The fit parameters 
are all dependent on x. c) Asymptotic value of packing den- 
sity poo vs x. The maximum occur at a; ma x — 0.9 while the 
jamming transition point at x g ~ 0.96. The dashed line rep- 
resents the exact equilibrium packing density. 




FIG. 2. Density profile p(h, x) at time t = 2 16 vs height 
h and vibration amplitude x. The bottom wall is located at 
h = 0. Here the particle- wall interaction is "attractive". 




FIG. 3. Density profile of a 1:1 binary mixture vs height 
h and time t. Top: unconstrained particles; bottom: con- 
strained particles. The vibration amplitude is x — 0.1. The 
system was prepared in a random homogeneous initial state. 
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